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Abstract 

The aim of this article is to justify mathematically, in the two-dimensional periodic setting, a general¬ 
ization of a two-phase model with pressure dependent viscosity first proposed by A. Lefebvre-Lepot 
and B. Maury in Adv. Math. Sci. Appl. (2011) to describe a one-dimensional system of aligned 
spheres interacting through lubrication forces. This model involves an adhesion potential apparent 
only on the congested domain, which keeps track of history of the flow. The solutions are constructed 
(through a singular limit) from a compressible Navier-Stokes system with viscosity and pressure both 
singular close to a maximal volume fraction. Interestingly, this study can be seen as the first mathe¬ 
matical connection between models of granular flows and models of suspension flows. As a by-product 
of this result, we also obtain global existence of weak solutions for a system of incompressible Navier- 
Stokes equations with pressure dependent viscosity, the adhesion potential playing a crucial role in 
this result. 

MSC2010 : 35Q30, 76N10, 76T20, 76T25. 
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1 Introduction 

In most of mathematical studies on the Navier-Stokes equations it is assumed that the viscosity 
is either constant or depends on the temperature and on the density in the compressible case. 
However, it is known that viscosities in real fluids, even incompressible ones, may vary not only 
with the temperature but also with the pressure. In his seminal paper [35] on fluid motion, 
Stokes already mentionned the possibility that the viscosity of a fluid may depend on the 
pressure. As explained by J. HRON, J. MAlek and K.R. RAJAGOPAL in [18] such dependence 
is for instance relevant for fluids at high pressures and flows involving lubricants. Another 
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example of pressure-dependent viscosities is provided by the theory of dense granular flows. 
The similarities shared by these types of flows with non-newtonian flows such as Bingham 
fluids, yield P. JOP, Y. FORTERRE and O. POULIQUEN to propose in (T9] a constitutive law 
based on a dimensionless number, the "inertial number" /. In their model, called "the p(I)- 
rheology", the volume fraction is linked to the inertial number T = <!>(/) and there exists 
a relation between the pressure P, the shear stress r and the shear rate D (u), u being the 
velocity of the fluid 

2d 

t = n(I)P with 1= -n-jr|D(u)|, 

( P / Qs )°- 5 

d and g s being respectively the particle diameter and the particle density. 

An important feature of granular flows is the existence of a maximal volume fraction 4>*, 
a constant approximative^ equal to 0.64 which corresponds to the random close packing. 
Taking account of such a congestion constraint leads to propose a model which can describe 
both the free/compressible regions where <h < <3?* and the congested/incompressible regions 
corresponding to = 4>*. The problem can thus be seen as a free boundary problem between 
the two subdomains. In dimension one, A. LEFEBVRE Lepot and B. MAURY in [22] proposed 
the following system which takes account of the previous constraints 

+ d x (<bu) = 0 
0 < < 4 >* 

< d t {$u) + d x {$u 2 ) + d x p = 0 (1) 

dtP a + d x {P a u ) = -p 
^ P a < 0, Pa (4?* - $) = 0 

The idea of the variable P a comes from [25] where it is seen as the adhesion potential of a 
single particle against a wall and measures in a certain sense smallness of the wall-particle 
distance. 

This article proposes to investigate a certain generalization of the previous system in the 
two-dimensional case 

' dt& + div ($u) = 0 (2a) 

0 < 4> < <f>* (2b) 

< dt($u) + div (<hu <g) u) + VII — VA — 2div (($ + II) D (it)) + r&\u\u = 0 (2c) 

<9tII + div (Uu) = (2d) 

,n>o, n($* - $) = o ( 2 e) 

Note that compared to we have three extra terms namely VII, — 2div((<h + n)D(u)) 
and r<&\u\u. The first two terms encode respectively the effect of some pressure law in the 
suspension model and the effect of the shear viscosity coming from the multi-dimensional 
setting whereas the last term, represents the friction. 

Following the ideas previously developped in HU and jSU, we approximate this system by 
a compressible Navier-Stokes system with singular (close to 4>*) pressure n £ . We also consider, 
and this is new compared to HD and ED, volume fraction dependent viscosities /r £ , A e singular 
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close to <h 


(3a) 


' d t $e + div ($ £ u e ) = 0 

< d t {<& £ u £ ) + div (<f* £ -u £ ® u £ ) + V7t £ (& £ ) + r® £ \u £ \u £ 
k -V(A £ (<f> £ )div u £ ) — 2div (/z e (<& e ) D (u e )) = 0 (3b) 

It is then expected that 7 t £ (<I> £ ) converges towards II, // e ($ e ) towards $ + II and A £ (<i? £ )div u £ 

towards A. As explained in [31 j . the singular pressure 7r e is not only useful for numerics, since 
it ensures automatically the constraint <I> £ < <!>*, but is also relevant from a physical point of 
view. It is indeed well-known in the kinetic theory of dense gases (see El) that the interaction 
between the molecules becomes strongly repulsive at very short distance. This effect comes 
essentially from an electrostatic force due to the fact that the electron clouds of different 
atoms or molecules cannot mix together. Several empirical formula have been proposed to 
describe this force (see for instance the general book [ITS] . or the famous paper [13] for a 
particular potential called the Carnahan-Starling potential), the common point of all of them 
is to consider singular potentials going to infinity faster than all the other forces involved in 
the model. Coming back to the theory of granular media, such repulsive pressures are also 
taken into account in the description of granular gases. In the gas regime, B. ANDREOTTI, 
Y. FORTERRE and O. POULIQUEN describe in their book [T] the kinetic theory that has been 
developed based on the principles of Boltzmann and Enskog. In particular some models involve 
the Carnahan-Starling potential. 

In the liquid regime of granular flows, one needs then to take account not only of the sin¬ 
gular pressure but also of a singular viscosity. Studying experimentally the case of suspensions 
(mixtures of fluid and grain in a dilute regime), one can define an effective viscosity of the mix¬ 
ture which is shown to vary with the volume fraction and which is expected to diverge close 
to the maximal volume fraction <f>* (see the books previously mentionned : [I] and |15|). As 
for the kinetic theory, only empirical laws are avalaible. From a mathematical point of view, 
to the author’s knowledge, there are few mathematical studies on fluid models with singular 
viscosities. However, one can cite the interesting paper [22 j. where A. LEFEBVRE-LEPOT 
and B. MAURY study a simple model in one space dimension of aligned spheres interacting 
through lubrication forces. From solutions of the discrete model, they construct a micro-macro 
operator and prove the weak convergence of the solutions towards global weak solutions of the 
continuous Stokes system (3>* = 1) 

, - a ‘ (tt$ s *“) = 

dt& + d x (Qu) = 0. 

To the author’s knowledge, this seems to be the first mathematical justification of the presence 
of a singular viscosity in a coupled system. 

For singular (close to <h*) viscosities and asymptotic two-phase description, it seems that 
nothing is known concerning mathematical justification. The problem has been envisaged by 
A. LEFEBVRE-LEPOT and B. MAURY in J22|. At the end of this paper they suggest that the 
singular system 

dt<& + d x (<&u) = 0 

d t ($u) + d x (<5>u 2 ) - d x ^ x ^ ^ d x v^J = <!>./' 
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could converge as e —>• 0 towards the hybrid system 


dt<& + d x (<bu) = 0 
d t {$u) + d x ($u 2 ) + d x p = / 
d t P a + d x (P a u) = -p 
k Pa<0 , $ < 1, P a (l - $) = 0 


previously presented. Nevertheless, the singular limit passage e 0 towards the hybrid 
Navier-Stokes system is not rigourously proven. Note that in the one-dimensional setting div 
and V are the same and also that they do not consider pressure in the momentum equation. 
This is the main difference between our mathematically justified asymptotic system (J2} and 
the proposed limit system (J4]). 

In this paper we want to take account of convection and of a singular pressure, a natural 
question is then to know if we have to impose a relationship between the singular viscosities 
and the singular pressure. An interesting remark for our study which can be found in [1] 
is basically the following: if one wants to describe within the same framework suspensions 
and immersed granular media (described by the jr(/)-rheology introduced before), one has 
to ensure the compatibilty of the two formulations by imposing the same divergence in the 
viscosity and the pressure close to <3?*. This is the approach followed in this paper where the 
shear viscosity p £ and the pressure ir £ increase exponentially close to <h* 


<J> 

= - ( exp 


1 —|—tt 


1 - <&/$* 


A e ($) = 2 £ ° 


<j > 2 


<f> 

<f>* | 1 - 

c[>* 


exp 


- 1 ) +*, 

£ l+a 

1 - $/$* 




x ‘ (4) = 7l^ 


7 


exp 




1 - <3?/<f>’ t 


- 1 


with <h* a fixed constant, a > 1 and 7 > 1 . 


The mathematical justification of the limit passage e —> 0 from a compressible model of 
type ® with a singular pressure ir £ towards a two-phase model of type ® has been the 
subject of two recent articles, EH and m respectively in the one-dimensional setting and in 
the three-dimensional setting with an additional heterogeneity in the congestion constraint. 
Nevertheless, these papers concern only constant viscosities and cannot therefore cover the 
case of dense suspensions for which we have seen that the viscosities depend on the volume 
fraction. To answer this question we need to carefully study the compatibility between the 
estimates derived from compressible Navier-Stokes equations with volume fraction-dependent 
viscosities and the singular limit passage e —>• 0. 

More precisely, considering "degenerate viscosities" (meaning that viscosity /u(<f>) vanishes 
on the vacuum, <f> = 0), one cannot deduce from the energy estimate a control on the gradient 
of the velocity contrary to the constant case. To deal with this difficulty, D. Bresch and 
B. Desjardins proposed in 0 , m a new entropy for the compressible Navier Stokes system 
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with degenerate viscosities. First for the shallow water viscosity /i(<3?) = /U 0< f>, A(<h) = 0 
(see also ®), then for more general viscosities /x(3>), A(<h) satisfying the algebraic relation 
A(3>) = 2 (//(<F)<f> — /i(<F)). The idea is to introduce the effective velocity w = u + 2V</?($) 
where <p is linked to the viscosity by the relation 






and to derive the energy associated to this velocity. In [19], D. BRESCH, B. DESJARDINS and 
D. Gerard-Varet proved the stability of the solutions for the compressible Navier-Stokes 
equations with additional terms such as drag terms or a singular (close to 0) pressure. With 
no extra terms, the stability result is given by a new estimate derived by A. MELLET and 
A. Vasseur in [28] . This estimate provides then the extra-integrability on \f$>u necessary to 
pass to the limit in the convective term of the momentum equation <8> u. 

In the shallow water case /r(3>) = <&, A(3>) = 0 with friction or cold pressure, D. BRESCH 
and B. DESJARDINS gave some hints in [ 8 ] to built a sequence of approximate solutions com¬ 
patible with the BD entropy and E. ZATORSKA in Wl [38] gave the complete proof of existence 
of weak solutions. Note that for the shallow water sytem with no drag terms nor cold pressure 
it is also possible to construct global weak solutions. The idea developped by A. Vasseur 
and C. Yu in [36] is to consider the system with drag terms and a quantum potential, then 
to construct a smooth multiplier allowing to get the Mellet-Vasseur estimate which does not 
depend on the drag. It is then possible to let the drag term go to 0 in the equations to recover 
weak solutions of the classical compressible Navier-Stokes system. 

For the general Navier-Stokes system with the algebraic relation A(<F) = 2(//(<!>)<I> — //(<!>)) 
construction of approximate solutions satisfying the energy and the BD entropy is not easy 
even with the hints given in [ 8 ]. In [10], D. BRESCH, B. DESJARDINS and E. ZATORSKA 
propose a new concept of global weak solutions called K-entropy solutions which is based 
on a generalization of the BD entropy. Considering the energy associated to the velocity 
w = u + 2 KV(p, where k € ( 0 , 1 ) is a fixed parameter, they derive the ^-entropy estimate. 
This notion of weak solution is weaker than the previous based on the energy and the BD 
entropy in the sense that a global weak solution of the compressible Navier-Stokes equations 
which satisfies the energy and the BD entropy is also a K-entropy solution for all 0 < k < 1. 
Yet, the opposite claim may not be true. 


In the framework of degenerate viscosities, to the knowledge of the author, the only work 
justifying the limit passage from the compressible system towards the two-phase model con¬ 
cerns the shallow water equations (/r( < F) = T, A( < I > ) = 0) with capillarity and a power law 
pressure a<l ? 7 with 7 —>• +00 (see [21]). Moreover the authors need to multiply the weak 
formulation of the momentum equation by to deal with the possible vacuum states $ = 0 . 

The objective of the present paper is to address first the global existence of weak solutions 
in the two-dimensionnal periodic setting to the suspension model |3]) with the singular pressure 
7 r £ and singular viscosities /r e ,A £ introduced before and an additional friction term necessary 
to ensure the limit passage in the convective term & £ u £ <8> u £ , as explained previously. In 
dimension 2 , imposing the same divergence on 717 and p, e we prove the global existence of 
K-entropy solutions satisfying the constraint 0 < < <F*. 
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The second part of the article consists of the justification of the singular limit e —>• 0 towards 
the two-phase system © modelling a granular media. Compared to the previous work with 
constant viscosities |31j . it is interesting to note that the singularity of the viscosity simplifies 
some compactness arguments and brings more regularity on the limit pressure. Indeed, the 
K-entropy gives then a control in dimension 2 of all the powers of p e and in particular, since 
we have chosen the same divergence on the pressure and the viscosity, this implies a control 
of the singular pressure 7T e with no need of additional estimate. 

All this study strongly relies on the uniform L°°(0, T; L P (I7)) controls, p £ [l,+oo) of the 
singular (close to T*) coefficients derived from the ft-entropy estimate due to the fact that the 
space dimension is equal to 2. 

As a corollary of our result, if II 0 > 0, taking as initial volume fraction <f>° = <f>* and 
approximating system © by © with an appropriate initial datum (u^, II 0 (#^), converging 
to (u°, II 0 , <!>*), we obtain weak solutions to the fully incompressible system 

' divu = 0 (5a) 

( dtu + u ■ Vu + V(II — A) — 2div ^ + l^j D (u)^ + r\u\u = 0 (5b) 

II > 0, <9*11 + div (ull) = — — (5c) 

2 

There have been few mathematical studies concerning incompressible flows with general 
pressure dependent viscosities. Most of the works, see for instance the interesting review paper 
m by J. MAlek and K.R. RAJAGOPAL or the article |12| . deal with a viscosity depending 
on both the pressure and the shear rate 

p = p(p,\B (u)\ 2 ) 

with an implicit relation between the Cauchy stress tensor and the shear rate D(u). It 
seems to be no global existence theory for purely pressure dependent viscosity. In [[32] . M. 
RENARDY confirms the physical relevance of a linear dependence of viscosity with respect to 
the pressure; indeed he proves that pressure driven parallel flow exists only if the viscosity 
is a linear function of the pressure. But M. RENARDY can establish (c./. [132]) existence and 
uniqueness of solutions only under a restriction on the velocity field : the eigenvalues of D (it) 
have to be strictly less than 

1 

linip^oo p'{p )' 

Later, F. GAZZOLA showed in m a local existence result without the previous restrictions, 
but for small data and assuming an exponential dependance of the viscosity with respect to 
the pressure . 

It seems then that there is no equivalent of our result in the literature on incompressible 
flows with pressure dependent viscosity. We obtain a global existence result of weak solutions 
with no restriction on the initial data nor an unrealistic relationship between the viscosity 
and the pressure. Although the ratio p/p tends in our study to 1 as p A +oo and not to 
Too as suggested in (23, our constraint is consistent with the arguments developped for the 
theory of granular flows in [Ij. Note the important role played by the adhesion potential II 
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in our mathematical results since it is this potential, via the equation (l5cl) . which provides 
the additional estimate necessary to prove the stability of the solutions of the incompressible 
system ©• 


2 The suspension and the two—phase granular systems 

As mentionned in the introduction our study restricts in dimension 2. To ensure the sufficient 
controls and the compactness of all singular quantities, in all the paper Q will be the periodic 
domain T 2 . We consider the two-phase granular system 


' + div ($> u) = 0 (6a) 

0 < $ < (6b) 

< c^II + div (n«) = — ^ (6c) 

dt($u) + div (<f>n <g) u) + VII — VA — 2div ((II + $) D ( u )) + r$\u\u = 0 (6d) 

_ $ II = <M > 0 (6e) 


with <3?* a positive constant which represents the maximal volume fraction and r > 0 a 
small coefficient which will be determined in the proof. We supplement the system by initial 
conditions 

®\t=o = $°, (*«)|t=o = m ° 

n° € l°°(P) n w 1 ’ 2 ^) (7) 

with 

0 < 3 >° < $* ( 8 ) 

= 0 a.e. on {x € ft ; d> 0 (x) = 0}, - € L l (Q) (9) 


Remark: As said above, the parameter r is small, precisely it is taken sufficiently small to 
ensure the compatibility of the system with the K-entropy inequality as it will be explained in 
Section 13.11 


2.1 A model for suspension flows based on singular Compressible Navier- 
Stokes equations 

We approximate the previous two-phase system by mean of a singular pertubation, we will 
call this perturbed system the "suspension model" 

' d t $ e + div (d> £ ii £ ) = 0, (10a) 

< d t {^ e u £ ) + div ($ e u £ <g> u e ) + Vvr £ ($ £ ) + r$ e \u £ \u £ 

—2div (/r £ (d> £ ) D (it e )) - V(A £ (<k £ )div u £ ) = 0 (10b) 
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The viscosities are defined by 


=Me($) + *= < 


1 | / c!+“ 

£ -*|ex P | 


- 1 


4>* / 


4> 

+ $ if — < 1 
4>* 


+ oo 


if 


i>i 

4>* 


with a > 1 and the algebraic relation introduced by BRESCH and DESJARDINS in [7J 

a £ ($) = 2 

The singular pressure is related to these viscosities and is defined by 


= ( Jr) a4($) = < 


$ / $ XT' 
£ \$* 

+ OO 


/ 


/ 


exp 


—. 1 —|—CL 


1 - 


<L 

4>* 


1 | if -5- < 1 

4>* 


4> 

if 7r>l 
4>* ~ 


with 


7 > 0. 


( 11 ) 


( 12 ) 


(13) 


(14) 


Remark: Formally we observe that a > 0 ensures the convergence to 0 as e —>• 0 of the 
singular terms /r*, A £ and 7r £ on the set {<L < 4>*} but we will see in the proof of Lemma [2] 
that we need a > 1 to guarantee the convergence of (1 — $ £ /4>*)7r £ ( < l> £ ) towards 0 and obtain 

m- 


The system fllOal) fllObl) is supplemented by initial data (<&[), m°) 


0 < < <f>* 


\m 


012 


4,0 


= 0 a.e. on {x € D ; = 0}, 


\m. 


012 


m 


e l\q) 






where e £ is such that 


€ L'i'A) 

e[ 1 4>2 ’ 

and where all the bounds are uniform with respect to e. 


(15) 

(16) 

(17) 

(18) 


Let us now introduce the notion of weak solution for system fllOal) - fjlObl) . 

Definition 1 (^-entropy solutions of ffTOl) ! Let T > 0 and k G (0,1), (<f> £ ,u £ ) is called a 
K-entropy solution to system fllOal) - fllObl) . under the initial conditions (fl5l) - fll8l) if it satisfies 


• rrieas{(f, x) : & e (t,x) > 1} = 0; 
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the mass equation in the weak sense 
r T r rT 


[ [ $ e dtZ- [ [ $ e « e -V£ = [ $°£(0) V£eP([o,r)xfi); 

Jo Jn Jo Jn Jn 


the momentum equation in the weak sense 


[ [ $> E u E -d t (- [ [ ($ E u E <g> u e ) : VC - / / 7r e ($ e )divC 

Jo Jn Jo Jn Jo Jn 


+ r [ f ® E \u E \u e • C + 2 [ [ ^ £ (<f> £ ) D (it e ) : VC 

Jo Jn 

f T 

+ 


to Jn 


f f A £ (<f> £ )divu £ divC = / m° • C(0) 

Jo Jn Jn 


for all (e (2?([0,T) x!l)) . 
the n-entropy inequality 

| u e + 2kV</? £ ($ £ )| 2 


sup 

tG[0,T] 

+ r 




+ k( 1 _ k) M^ |+ ^ ($e 


sup [ n E {<S> E ) + r [ [ $ £ |u £ | 3 

e[o,t] Jn Jo Jn 

+ k/ T / M,(-f e )|A(» e )| 2 + 2K [ T f 
Jo Jn Jo Jn 


|V$ £ 


+ (!-«)/ [ M$e)|D(u £ )| 2 + f (/i'($ e )$ £ - /x e ($ e )) |divu £ | 

jo uii J n 

<C(r) 

where ip E is such that 


til*) = 


til*) 


$ 


and C(r) is a constant which depends only on r and on the initial data ($J,ra®) 


(19) 


( 20 ) 


( 21 ) 


Remark: D (it) and A (it) are respectively the symmetric and the antisymmetric parts of the 
gradient defined by 

^ , . Vit + V*u . , . Vit — V f u 

D («) =-j-’ A(u) = - 2 -' 

Remark: To prove the stability of the solutions the integrals of the diffusion terms 

2 / / Me(3> £ ) D (n e ) : VC + / / A £ ($ £ )div (u £ )divC 

Jo Jn Jo Jn 

should be understood in the following sense 


diHel*e) 


Ve{*s 


! / / /*,(*„) D (u«) : V C = - f f vT<4 ( + PP4J 

jo Jn Jo Jn 


f / 

Jo Jn 


i (djhe($e 




V$~e 

V + 


( 22 ) 
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and 



A e (<f* e )div u £ divC 



X / C ^ £ J , VU £ div C- 


(23) 


Remark: Compared to the work of BRESCH, DESJARDINS, ZATORSKA [TO], we have the 
additional integral 

r sup / He($e) 
te[o,T] Jo. 

which allows us to deduce the control of /i e ($ £ ) in L°°(0, T; L p (fl)) for p € [l,+oo) if we 
combine it with the control of V/r e (<h £ ) in L°°(0, T; L 9 (U)) for all q € [1,2). Note that this 
control strongly relies on the presence of the drag term in the equations. 


Remark: It would be also possible in the analysis to consider the more simpler case of 
pressure and viscosities singular close to <h* as power laws, namely 


7T e ($) 


e<f* 



P 


as proposed in [22] . keeping the relationship between the coefficients p £ , ir £ and X e . The 
existence of weak solutions when e is fixed works exactly in the same way. Modifying slightly 
the arguments, one can also prove the limit passage £->0 towards the same hybrid model 
satisfying (<h* — <h)II = 0. 


2.2 Main results 

Under the conditions previously stated, we are able to build global weak solutions of the 
system fllOal) (jlObj) . 

Theorem 1 (Existence for the suspension model) Let T > 0, £ > 0 and (<k[! ,771°) an 
initial data satisfying (USD-dUD. There exists r > 0, which depends on T, such that there exists 
a n-entropy solution (<h e ,u e ) to the suspension model fllOal) - fllObl) in the sense of Definition 

m 


Thanks to the previous existence result we can address now the question of the singular limit 
passage e —>• 0 towards the two-phase system. 

Theorem 2 (Existence for the two phase system) Let.T > 0 7 (<h 0 ,m°,n 0 ) and ($g, m®) 
satisfy respectively (0-® and (1151) — (|18l) . We assume that —> l* 0 in L P (D) for all p € 

[1, Too), mfi e j yTf>0 —» in L 2 (fl) as e —>• 0. Let r such that we can apply Theorem 

[7] Then there exists a subsequence (<I? e , u e , ir e (<& e ), A £ (<f> e )div u £ ) converging to (T>, II, A) 
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solution of 


f [ mi - [ T [ $u ■ v? = [ $°e(o) 

Jo Jn Jo Jn Jn 

0 < $ < 4>* a.e. in (0 , T) x fl 

[ T [ ud t c- l T [ n«-ve + <4^)= / n°e(o) 


/ 0 JS! 
r T 


1 0 J£2 

rT 


in 

rT 


(24a) 

(24b) 

(24c) 


— f f &u-dt(— f f (4>u<g) u) : V£ — / / Ildiv^ + r / / • C 

jo jo in jo Jn jo in 

+ (A, div C) + 2( (n + $) D (u), VC ) = f m° ■ C(0) (24d) 

Jn 

4*11 = 4**11 > 0 a.e. in (0,T)xSl (24e) 


for all £ £ 2?([0, T) x fl), £ € (D([0,T) x II)) 2 and where the terms (— , £) and (A, div£) /iaue 
to be understood as the duality pairing between the distribution A and the test functions £, 
div£. Moreover, the limit has the following regularity 


$ € C([0, T|; ^(fl)) n I°°(0, T; ^ ll2 (fl)), for all l<p<+oo, 

n £ L°°(0, T; W 1,2 (M)), 

VTTT® d (u) € a 2 ((o, t) x n), 

V&u € L°°(0, T ; L 2 (J4)). 

A £ W _1,oo (0, T; L P (I2)) + L°°(0, T; W~ 1,q {M)) for all p £ [1, +oo), q £ [1, 2). 


Remark: We observe that we get much more regularity on the limit pressure II than in the 
constant viscosities case |31j . As we will see in the proof, this is a consequence of the K-entropy 
and the relationship satisfied by fi e and tt £ . In particular this regularity gives a sense to the 
product 


(IID(«), VC) 



(25) 


The difficulty in the proof of Theorem [21 compared to the case £ > 0 relies on the fact that 
at the limit e = 0 we do no have meas {( t , x ) : 4>(t, x) = 4>*} = 0. We then need to carefully 
study the control that we have on the singular coefficients taking into account the possible 
convergence of 4> e towards 4>*. 
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Global existence of weak solutions to an incompressible Navier-Stokes system with 
pressure dependent viscosity 


As in the constant viscosities case studied by P. L. LIONS and N. MASMOUDI in | l24] . we 
prove in Section [5] the compatibility on the limit system between the constraint f!24bl) and the 
divergence free condition 

divu = 0 a.e. in {<£ = <!>*}. (26) 


If initially the two-phase system is entirely congested, meaning that <b 0 = <3?*, II 0 > 0 and 
divu 0 = 0, then, considering the approximated singular system fHoa-fflObl) with initially 
n ( $*\ 

<j>u = <|>* / i _ £ a -__ 1, the previous theorem will give us the existence of global weak solutions 
for the incompressible system with pressure dependent viscosity. 


Theorem 3 (Existence for the incompressible system) Let T > 0, (u°,n°) such that 
u° G L 2 (Q), divu 0 = 0 and n° > 0, with 

n° g L°°(n) n w li2 (fi). 


There exists a global weak solution to the pressure dependent incompressible system for all 

£gp([o,t) xfi), Ce mo,T)xn)) 2 


div u = 0 


(27a) 

(27b) 


satisfying 


- [ T [ mz- f T f nu-v£ + (^,£>= f n°e(0) 

Jo Jn Jo Jn A Jn 

f [ u-dt c- [ T [ (u®u):VC- [ T [ -j^divC + r T [ \u\u ■ C 
Jo Jn Jo Jn Jo Jn ^ Jo Jn 

+(-^,divC)+ 2 ((II/<r + 1 ) D(u),VC) = J u°-C(0), (27c) 
n>0, n€L oo (0,T;W rl ’ 2 (f2)), 

A G W 1 ’ 00 ^, T; I?(Q)) + L°°(0, T; W 1 ’^)), Vp G [1, +oo), q G [1,2), 
u G L°°(0, T; L 2 (Q)) D L 2 ( 0, T; W 1)2 (fi)). 


Remark: We can recover the regularity of the potential II directly from the system (I27al) 
(I27cl) using the BD-entropy (see for instance [6]). To simplify the explanation, we drop the 
drag term of the momentum equation. Taking the gradient of 

d t U + div (Uu) = — ^ 


and dividing by <!>*, we have 


9tV (^ + 1 ) +diV ( W0V (^ + 1 )) +diV 




VA 

2d>* 


= 0 
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Then, introducing the effective velocity w = u + 2V + 1 ), w satisfies, 

vn // n \ \ VA VI 

dtw + div [u <g> w) + —— 2div ( ( ^7 + 1 ) A (u) ) + -J- - -±— = 0. 

Finally, multiplying this last equation by w and integrating, since u is divergence free, we get 

2 


\ISa ?+ L 


vn 




+ 2 l n li^ + 1 ) |A ^ _0 


which ensures that w is in L°°(0, T; VF 1 , 2 (fl)) if initially u° £ L 2 (fl) and n° £ VF 1 , 2 (fl). 


2.3 Sketch of the proof of Theorem |T] 

The proof of Theorem Q] is not a direct consequence of the theory of the K-entropy developped 
in m since the pressure and the viscosities are singular close to <3?*. To deal with this 
difficulty we first add a parameter 5 in order to truncate these singular terms. Then we add 

<J) 2 

an artificial pressure $Vp(<3?) = i9V —,'d > 0 in order to control the gradient of the density. 
The approximate system reads as 

' dt®^& + div ($^5U#,6) = 0 (28a) 

+ div ($^, 51^,5 (8) u& >s ) + + Vtt^'F,^) + ^#, 51^,6 

— 2 div (He,s(®0,s) D ( u #,s)) - v (A e , 5 )div (u^g)) = 0 (28b) 


with 


tt £ ,5(^) = < 


T / <I> 


£ V <f> 


exp 


$ / <f> \ 7 / 

7 11 7 * * 1 


exp 



<j> 

- 1 I if — < 1 - <5 

<j>* 


$ 

if — > 1 - 6 


(29) 


VeA®) = < 


$ 




exp 


V 


1 - 


$ 

cj>* 


— 1 | + <J> if — <1 — 5 

cj>* 


l+a' 

exp | —) - 1 ) + $ 


$ 

if — > 1 - 5 

cj>* 


(30) 


£ v " V s 

If we denote (fJ> e s)+ the right-derivative of /j, e g, \ £t $ is related to // e § via the algebraic condition 

a e ,«($) = 2 (0m)'+(*)* - am( $ )) 

( 


= 


2e a - 


& 


<F* 1 - 


<f> 

<I>* 


■ exp 


£ 


1-f - ^ 


$ 

(|>* 


$ 

if — <1-5 
(|>* 


<h 

if — >1-6 
< 5 * 


(31) 
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The proof of Theorem [T| will consist of two limit passages, first 5 —0 then i? —>• 0. 

Remark : One could try to extend the procedure of VASSEUR and Yu in [36] and derive a 
control on <3?|rt| 2 log(l + |w| 2 ) but it seems not possible to get a uniform estimate with respect to 
e or even 5. This is the reason why we need to keep the the turbulent term in the momentum 
equation (|28b|> . 

Organization of the paper 

The rest of this article is coarsely divided in three parts. The first part concerns the proof of 
Theorem |T] namely the existence of weak solutions for what we call the "suspension model" 
fllOal) fllObl) with singular viscosities p e , A e and singular pressure n E . The main part, corre¬ 
sponding to Theorem [2] consists of passing from solutions of this suspension model towards 
solutions of the two-phase system of granular type f!24al) - fl24el) . Finally using this result we 
approximate the incompressible model (I27al) ~ (|27cj) by an appropriate suspension system and 
prove therefore the existence of global weak solutions for (I27al) - fl27cl) as stated in Theorem [3] 


3 Existence of solutions to the suspension model 

3.1 Global existence of K-entropy solutions when e,5 are fixed 

We first need to prove the approximate system containing all the parameters e, 5, admits 
global weak solutions. We recall in the following definition the notion of K-entropy solutions 
for system fl28al) - (|28bl) . 

Definition 2 (^-entropy solutions for fl28l) l Let T > 0, k € (0,1), is called a 

if it satisfies 


(32) 


n-entropy solution to system (128a[) 

• the mass equation in the weak sense 


- [ T [ $6dtZ - r [ $5U5 • V£ = [ $°£(0) € V([0 ,T) x 0) 

Jo Jn Jo Jn Jn 

the momentum equation in the weak sense, € (D([0, T) x fl)) 2 

^sus-dtC- [ [ (®5Us ®us) : VC - d [ [ p($s )divC 

Jo Jn Jo Jn 

n e ,6{$s)div( + r / 4> 5 |u 5 |?x 5 • C 

Jo Jn 

[ [ He,, 5(^a)D(u 5 ) : VC+ / / A £i5 ($ 5 )div(n 5 )divC 

Jo Jn Jo Jn 

= [ m° e ■ C(0) 

Jn 


to Jn 

rT 


If 

rT 


(33) 
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• the K-entropy inequality 


sup / $5 
te[o,r] Jfi 


\us +2KVip £t s(^s)\ 2 , ,_J2V< / 9 £i ,5($5)P 

--- | /^( 1 Kj) --- 


+ sup [ f^^ + $ 5 e £!5 ($ 5 ) + r// £ia ($ (5 )') 

te[o,T] Jn V 1 J 

+ «jf/ n Mm(<MI A Ml 2 + 2 kJ q T f Q <*(*«) ^ |V4> 5 | 5 

+ (1 -k)[ f H e j{$6)\~D (u S )\ 2 + [ -/^ £ , 5 (^« 5 )) |divu, 5 | : 

Jo LJn Jo 

[ T [ <*> £ M 3 < C(r) 

Jo Jn 


+ r 


where tp £ 5 is such that 


<M*) = —*—■ 


(34) 


Remark: C(r) is a constant which depends only on r and on the initial data (<h £ ,m°). As 
we will show later on, this constant will derive from a non-linear Gronwall inequality. 

In [10], D. BRESCH, B. DESJARDINS and E. Zatorska base their construction of approxi¬ 
mate solutions on an augmented approximate scheme satisfied by ($>,w = u + 2 kV(/ j( < h),u = 


2V</?(<3?)). In our framework this augmented system writes as 

' d t ®s + div (4 >sws) - 2nAp, E} s{$5) = 0 (35a) 

dt(®sws) + div ($sus ® ws) - 2(1 - K)div (/i £ij (fj)Vwj) - 2Kdiv (p £! s($s) A (w 8 )) 
+i9Vp(4> 5 ) + V7 t £i 5 (<I\5) + 4(1 - K)Kdiv (/^(-h^vV^^)) 

-V((A £ , 5 ($ 5 ) - 2 k(h' £ s(®s)$6 - ^ e ,j( < f , 5 )))div u 5 ) + r^ 5 \u 8 \u s = 0 (35b) 

< 

^(4> 5 V^ £;5 (^> 5 )) + div ($ s u s ® Vip £j s(^s)) - 2Kdiv (p £ A®5)^ 2 <Pe,s(^s)) 

+div (/Lt £i5 ($ 5 )V 4 u; 5 ) + V((/4,5($<0$<S ~ fi e ,s{®s ))div u s ) = 0 (35c) 

, ws = us+ 2KVtp £t s(®8) (35d) 


Justification of the K-entropy inequality : Following the steps of D. BRESCH, B. DESJARDINS 
and E. ZATORSKA |10j , we multiply ()35bl) by w 8 = u 8 + 2kV<£> £ 5 (^, 5 ) and we combine it with 
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the equation (I35cl) tested by 4/c(l — K)V<y9 £j ,5(<I>5), we get 


$<5 


I u s + 2 kV<^ £ 5(^)1 


+ 


+ 


in l 


+ $s e £ , s ($ s ) 


+ k(1 — k) 

(t)+r [ [ $5KI 
Jo Jn 


12V </7 e ,<5 (^< 5 ) | 


2\ 1 


(t) 


2 »/‘/^,)W^ + 2 K /7 ^ ( ^ ( * >) |V« < |» 

jo Jn Jo Jn **5 

+ 2(k — 1) f f ^ e ,5($«5)|D(u 5 )| 2 + 2(k- 1) f [ (n' eS (^s)^5 - VeA®s))(dwu s )' 
Jo Jn Jo Jn 

n (+ 2kV^(3> £ °)| 2 , |2V ¥M (4> £ °)|‘ 


+ k(1 — k)- 


(36) 


+ 


^Ml! + $o ($ o) 


-2 nr [ f $ s \u8\u5-'V l Pe,5(®6) 
Jo Jn 


In addition, the renormalized continuity equation writes 

dt(^e,s(^s)) +div(/r £i5 (4> 5 )u (5 ) = -A £i5 (4> 5 )divu 5 
which gives, after multiplication by r and integration in time and space, 


' / /i e ,i($<s)(t) + r [ [ 
Jn Jo Jn 


^e,s(^s) A . n 

—^-r-divu 5 = 0. 

n 2 


(37) 


Adding then (1371) and fl36l) 

I us + 2kV(^ £) 5(4> 5 )| 2 




+ k( 1 — k) 


+ 


+ 


<J)r 

tf-y + ^<5e £i5 (4> 5 ) + r/i £;5 (4> 5 ) 


12V (*^*6) I 
2 
rt 


2 \ 1 


(t) 


(t) + r 




0 Jn 


2k n + 2k ff 

Jo Jn Jo Jn ^5 

+ 2(k — 1) /" [ am($ 5)| D (u 5 )| 2 + 2 (k — 1) [ [ (fj? s ($s)®8 ~ AM(3>5))(divu 5 )" 
Jo Jn Jo Jn 

= f \*o(h°e/® 0 e+^\7<PeA$e)\ 2 , j2V¥> Ei a($°)| 5 

Jn 


+ 


i£lL + $“e M ($«) 


2 


+ k(1 — k)- 


— r 


ff 


K.si^s) a . 

—-div u 5 

n 2 


- 2 nr [ f d> 5 |u 5 |u 5 • V^a^a) (38) 

Jo Jn 

We now have to control the last two integrals of the right-hand side in the previous relation: 

J\ = 7 T / [ A £i a($a)divua, J 2 = 2nr f f $aKK • V<p £j a($a)- 
2 Jo Jn Jo Jn 
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Control of J\ = | j[J f (} Ag^^jdiv U 5 . Unfortunatly it not possible to control uniformly with 
respect to all the parameters directly by the left-hand side of (1381) . The idea is to apply in 
a certain sense the operator (—A) _1 div to the momentum equation where A ” 1 denotes the 
inverse operator of the Laplace operator. For each function / such that Jq f = 0 we denote 
g = (—A) -1 / the unique periodic function such that —A g = / and f^g = 0. Therefore we 
can obtain the equality 


rt r r f , . r 


r [ [ " £ ’°^ 0 / di vu 5 = - I (-A) 1 div($ 5 w 5 )(t) - - / (—A) i div(m^) 


r 

2 ./o Jn 

t 


A 1 d i dj(^sus u s) 


r 

+ 2 * 


+ r 


[ [ pm 

Jo Jn 

f [ A -'did. 
Jo Jn 



+ 9 / / ^, 5 (^ 5 ) 

z Jo Jn 


He,s($6) 


di u 5 + d j u S 


„2 r t 


+ 


f [ (-A) Miv^l^ln^) 

Jo Jn 


Ik 

k =1 


(39) 


In order to fully justify the previous equation, it suffices to take in the weak formulation (1331) 
the test function 


CiWVA- 1 ^)), a ec°°((o,T)), C2 € c°°(n), [ Ux) = o. 

Jn 

We refer to [16] Section 2.2.6 for a similar computation and to Section 10.16 for properties of 
the singular operators involved in the previous equation. 

Before studying each integral I ^ of the previous equation let us explain how the norms of 
^£,5(^5) are treated. By the Gagliardo-Nirenberg inequality, see for instance |29] (Theorem 
p.12 with j = 0 and m = 1 and Remark 5.), we have 


< C (ll/ i e,5(^ > 5)|lLi(r2)ll V / i £,<5( < ^<5)lli g -(% + 11 Me, 5 (^ 15 ) Hi 1 (Q) 

for all q E [1, 2), 


qe 1 


and where 


We obtain then 


2 q 

2-9 

1 c . l~g 

q 29 /( 2 - 9 )' 


IIM£,5(^5)IIl9(o) < C'l||MM(^ , 5)llii(n) + C 2 II— ^=^||L2(Q)||A/^'|| L 2 9 7(2- q -)( f2 ) 


< C , 1 ||^($5)||“ 1 i (n) + 


ai , 


^ \/d\5 


L 2 (n) + C , 2^||M£,5( < ^5)||L9/(2-q)(n) 
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If we consider now q such that 


q€ 


q 2q 

2 -q' 2 -q 


and v > 0 small enough the last term can be absorbed by the left-hand side and replacing 
He,s(^s) by THe,s{^s) we deduce that 


lkMe,5(^5)||L9(n) < C 


\\THeAm a L \ m 


+ 


^He,s($s) ||a 2 


(40) 


with «i and ct 2 two positive constants depending only on q. We come back now to (1391) and 
detail the control of each integral . Thanks to the regularity properties and conservation of 
mass we have for I\ (and the same for I 2 ) 


\h\<r |( —A) 1 dlV (^ 5 Us)\ < Cr\\y^\\ L oo L 2\\y/^ u 5\\L^L2 ^CrWy^usW^Li 

Jn 

which can be absorbed by the left-hand side provided that r is small enough. 

Concerning the convective term in I 3 , we use the estimate due to the drag and the control of 
this control coming from the splitting 

ll$5 +1 ||L 1 (fi) = ll^ > 5 +1 ( 1 {$«<3>*/ 2 } + 1 {<!>6>$*/2})IU 1 (n) 

< C + C[|^>5e £i 5(<h < 5)|| i i(n) 


Hence, 


\h\ < r 


fS 

Jo Jn 


A 1 did j ($ 8 Usu 3 s ) < Cr\\A l d i d j {^ 8 \u 5 \ 2 )\\ L i L 3 (- 1 +i)/( 3 + 2 1 ) 

<Cr||$J / 3 ^ / 3 |^| 2 || Lli 3 ( 7 + i)/( 3 +27 ) 

<^ 1/3 ||^ 3 || L3L3(7+1) ||r 2 /3$f |^| 2 lli3/2L3/2 

<C'r 1 / 3 ||^||^ T+ 1 ||r 1 /3$y3 U 5 | |2 3L3 
< Cr 1 / 3 ^,5^6)11^ 1 Ik ' ^<5 U S\\ 2 L 3 L 3 

Concerning the pressure terms, they have a positve sign and we do not need to control them. 
The integral of the drag can be controlled as follows 

|Cr| < r 2 [ [ (-A) _ 1 div($ 5 |n 5 |?X5) 

Jo Jn 

< Cr||r 1 / 3 ^ /3 r 2 / 3 $ 5 /3 |n (5 | 2 || i i i9 for all q > 1 


< Cr\\r 1 / 3 Ve,s(® 6 ) II L 3 L 3 q /( 3 - 2 q ) \\r 2 / 3 <f> 2 / 3 \u s \- 1| L s/ 2L s/ 2 

r-V 3 (J>y 3 


< Cr ( ||r^(^)HS L , + II ' !) ||g L , 


5 u S\\l 3 l3- 


Finally, the last integral Iq, coming from the viscosity He s(^s) D (u$), is the most difficult 

rt r r /diy? s + dju\' 


h = r 


f / A ~‘ 


didj 


He,s($s) 
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We see then that it is sufficient to control the gradient part, namely 

r J In A ~ ldidj 

the remaining term can be then treated exactly in the same way. Writing that 

^e,s(^s)9iU 3 g = di(n et s(®8)Us) - di(fj, et s(^s))u 3 s 

we observe that the first term combined with the operator A ~ 1 didj and the integration over 
Q will give us 0. We get 

r J J A — 1 cl ,;dj (ji EiS (^ s )diu{^ 

< Cr\\\/fi £> s(^s)u5\\ L ^Li 

< Cr 1_1/3_1/6 ||^^^|| i 3 L 2||r 1/3 d>y 3 M 5 || L 3 i 3||r 1/6 d>J /6 || L 3 L 6 9 /(6-5 9 ) 

< Cr 1 / 2 \\ V ^^^ 5 \ \L3L2\\r 1/3 ^] /3 U S \\ L 3L4 r ^sA^s)\\ 1 /3 L q/{6- 5 q) 
for all q € (1, 6/5) and fl40l) yields 

r J J A~ 1 d i dj ([i Sl s($s)diu{j 

< Cr 1 ' 2 (||r W (* J )|© 11 + || V '‘^5 >i! ||g t ,) ||r 1/3 ®J /3 «j|| L i tI . 

This concludes the control of the integral Jq/q A £! s(&s)div us in (13811 . 

Control of J 2 = 2 nr • V<p £i s(^s)' It remains then in fl38l) to control the integral 

coming from the drag 


I = 


2 nr / / |it5|w,5 • 

Jo Jn 


= - 2 Kr / / /i e ,5(^5)|^5|div (u s ) - 2 ur 

Jo J n 


[ [ Ve,6($s) 

Jo Jn 


u 




Splitting \ 7 u$ between its symmetric and its skew-symmetric part, we get 

\I\ < c(Q)nr [ [ |m £ ,5( $ 5)I|V«,5|KI 

Jo Jn 

<c(Cl)nr ( [ \/j, £> s($s)WD (us)\\us\ + c(Q)nr [ [ |/r £)5 ($ 5 )|| A (u 5 )||« 5 | 

Jo Jn Jo Jn 

y 2(1 — K) $ J 


+ c(CL)yfKr 2 ^ ||- ^ e 'y 3 - ^( U s)\\L 2 L 2 \\r l/i ® l s*Us\\ L 3 L -3 

<C' 1 r 2 / 3 ||V2(l-«)^MM(^)D(u 5 )||| 2L2 + C 2 r 2 / 3 ||^^y$yA(u 5 )||| 2 ^ 

+ C 3 r 2 l 2 \\r'l^'l\ s \\l, L , + Cr ^^ Wl... 


(41) 
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For r small enough we can ensure that 


max (Ci, C 2 , C 3 ) r 2 / 3 < ^ 


and absorb the first three terms of I by the left-hand side of (1381) . Splitting the last term into 
two parts 


A Pe^i^s) n _ ^ I, y/^e,s(^s) 1 


<F 


1/3 


L6L6 < II-^75-l{$i< 4 -*/ 2 }lli 6 i 6 + II- -^T /3 - 1 {'I > i>$*/2}ll^ 6 i 6 


-777-U 


$ 


we have on one hand 


y/»e,8(*s) 1 


<F 


1/3 - L { 4 ’ 5 <$*/ 2 } - 


1/6 


exp 


--1 —(—<11 


- 1+1 


L {^<®*/ 2 } 


< c$y u i{$ 4 <$*/ 2 } 


which is bounded in L°°((0, T) x fi) and on the other hand, 


APeA®s) 


Using again fldOl) 

llv'^,s(*j)i(» J >f/2>llfL‘ £ c (|I/m(<MII2!>li + ll v ^* ,) l 

Finally, coming back to (1381) . there exists a > 1 and a r > 0 such that 
1*45 + 2kV(/? £i(5 ($5)| 2 |2^ £i5 ($ 5 )| 2x 1 


Q2 

L 12 L 2 


+5 


it) 


+ 


+ k (1 — k)- 

[ [^«5e £i5 (+ 5 ) + r/j, ejS (®s)](t) + 7 ; [ [ ('dpi®s) + n e>s (®s)) + 
J ^ J 0 J al 


®5\us\ c 


0 1 !) 


+ K 


^s(^K,s(®s) 

®8 


|V ^| 5 


+ (k — 


[ [ Pe,8i®5)\A{us)\ 2 + 2 k [ [ 

J 0 Jn J 0 Jn 

1) [ [ Pe,5(®6)\V(u 5 )\ 2 + (k- 1) [ [ (v' £}S (®5)®8 - PeA®8))(dwu S f 

Jo Jn Jo Jn 


< 

Jn 

+ Cr ar 

x 




|mg/+g + 2KV^(d>g)| 2 + k(i _ ;:) |2Vy e , a ($°)[ 2 < 


+ 


In L 


l 0\2 




+ ^e £ , 5 (d> £ u ) 


0 \Jn 


®s 


\u s + 2n\/cp £t s{®s)\ 2 , M j 2 V<p e ,*($a )| 2 1 , ^ 

---h k(\ - k) --- + $5e £i< 5(^>,5) + r/j, £ ,s{®s) 


For r small enough we can use a nonlinear generalization of the Gronwall Lemma on the 
interval [0, T], see for instance [4] or 0 (L emma II.4.12 p.90, with f(y) = r ar y a ). 

Hence we have closed the kappa-e ntropy inequality which allows us to use the existence result 
of D. Bresch, B. Desjardins and E. Zatorska [TU] . 
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Proposition 1 Let T > 0, e, 5 be fixed, there exists r > 0, depending only on T, such that 
there exists a global K-entropy solution to system (1281) . In particular we have the following 
regularities 


V$s\ws\ € L°°(0,T;L 2 (Q)), 
eL°°(0,T-L\n)), 


V^\V^ £ ,S^s)\ €L°°(0,T;L 2 (P)) 
$S e L°°{0,T-,L 2 (n)) 

r$s\us\ 3 e L\(0,T) xfi) 


yJn e ,s($8)Vu s € L 2 ((0,T) x 0), ^/A e ,i($a)divu 5 € L 2 ((0,T) x fi) 

V<f> 5 € L 2 (( 0 ,T) x 51). 


N 




(42) 

(43) 

(44) 

(45) 

(46) 


3.2 Proof of Theorem[l} Existence of weak solutions for suspension model (TTT71) 

We aim here at proving Theorem [Q by letting the parameters 5 and D go to 0. 

More precisely, we try in this section to derive the uniform controls dealing with viscosities 
and pressures which become singular as 6 —>• 0 and using the friction term to ensure the 
compactness of the approximate solutions. Passing to the limit 6 —> 0 in the equations, we 
prove that the limit volume fraction satisfies the maximal volume fraction constraint 

0 < & £ ^(t,x) < d>* a.e. on (0, T) x Ll 

Finally we perform the limit passage D —>• 0 which means that we eliminate the artificial 
<t > 2 

pressure i?V—. This step does not present additional difficulty and will be briefly explained 
in the final remark. 


Estimates 

We recall that the singular terms write as 


Ke,s(®8) = < 


$5 /£p 7 

eld)* 


l ( 1 + \ 

pl+a 

exp ' 


V 1 d >*) 

£ 


$8 


- 1 


$6 f*sY / /cl+a 


£ V**J l6XP 


- 1 


if ^ < 1 - 5 

<j>* 


if 


cj>* 


>1-5 


he,s(^s) = < 


$1 

£ 


exp 


1 —l - d 


1 - 


d>* 


- 1 I + $s if S < 1 - <5 


1 +a' 

exp I —— ) - 1 ) + d > 5 if 


d>5 / t £ 


(J>* 


>1-5 
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K,s(®s) = 2 ((^£,5)+( < ^<5) ( I > 5 - Me,a($5)) 


= 


2 e a 


*6 


-. 1 —|—ct 


$* 11 - — 
4>* 


■ exp 


1 - 


$1 

cf>* 


if 

<j>* 


if 


— > 1-6 

(J)* ~ 


Control of &s- Thanks to the K-entropy inequality and to the bound /if s ($,5) > 1 we ensure 
that 

V \Z$g is bounded in L°°(0, T ; L 2 (fl)) 
and then that <£,5 is bounded in L°°(0, T; L p (fl)) for all p G [1, +00). Moreover 


II^V //U ^ (<&<s)|v ^ lllL2i2 - c 

which means that is bounded in L 2 ((0, T) x fi). 


Control of &s u 5 - For the momentum we directly have &s u 5 = V&sV'&sUg bounded in 
L°°(0, T; L q (f2)) for all q G [1,2) and bounded in L 2 (0, T; L 1 (n)) by writing that 

V($ a ua) = + 2y/W s u s ® 

In addition, the drag contribution the K-entropy inequality provides 

& 5 U 5 = bounded in T 3 (0, T; L q {Ct)) V</G[l,3). 


Controls of the viscosities. Thanks to the K-entropy we control uniformly -\/<$>gVip et g(<&g) in 
L°°(0, T; L 2 {CL)). If we set such that 

= v^U^) = ^4=^ ( 47 ) 

V ^<5 

we deduce that VI4 5(^5) is bounded in L°°(0, T ; L 2 (Q)), coming back to the viscosity Me 5(^5); 
using the bound L°°(0, T; L P (H)), p < +00, on we obtain 

V/i e ,i($ 5 ) = bounded in L°°(0, T; L*(fi)), Vq G [1,2). (48) 

Since we control p, e g($>g) uniformly in L°°(0, T; L 1 (II)) we deduce that 

he,s{^s) bounded in L°°(0, T; L p (fl)), VpG[l,+oo). (49) 

Me,5(^5) D (ug) bounded in L 2 (0, T; L q (£l)), V(/G[l,2). (50) 

We next bound the other viscosity coefficient A e ,s(3>,s) by comparison with 

--5 

,L°°{ 0 ,T-LP(n)) p<+ 00 

-6 


Me,5(^5) = < 



if — < 1 

4>* 


if 


h>, 

<J>* — 
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Indeed, the ratio for < I ) 5 / ( h* < 1 — 6 reads as 

_l+a 


^£,$($< 5 ) 


2e a <hc ( 


* exp \T^J¥ 


{ nl 5 ($ 5 )) 2 

Using the fact that 




-.l-j-fl 


exp 




1 - ^Z^* 


r l-|-<2 


we get then 


= exp | 

VI - $<5Z$ 

> exp | 

/ £ l+a 

A-^s/^ 


£ 2 ( 1 +a) 


6 I ex P 


1 — exp — 


1 - 3>s/®* 


- 1 


2 ' 




i - Qg/®* 
Vl-exp (- £ ‘+“)) 


A e , 5 (^, 5 ) 

<C{£) 


exp 


r l -|~CL 


1 - Qg/®* 


because X 1 —> X 2 exp(-X) is a bounded function on (0,+ 00 ). Since e is fixed at this stage, 
we can conclude to the control of X £ g 


A,s(®s) bounded in L°°(0,T; L P (Q)), Vp € [1, + 00 ) 


(51) 


and in addition 


X / £ ' S ^l\ = Xc ' 6 (l S \ 2 l J ‘ s ’ S ^^ 3/2 boundedin U°°(0,r;L p (U)), Vp € [l,+oo). (52) 

vVmI^) He A®*) 2 


To pass to the limit in the diffusion terms written as in 

quantity LAAA. L e t g ( 0 , 1 ), on the set {^Z^* ^ ^} 
V$5 


, we need also to control the 


He A® s) K HeA^s) 

Ah - Ab/$* 

is bounded in L°°(0, T ; L P (Q)), p < 00 . On the other set {d^Z^* < 3>}, <£,5 is far from <h* and 
then 


HeA^s) , Ah ( / £ 

< -I exp 


1 ~|-CL 


-1 ) +A$~6<CA$~S 


Ah - e V VI -*s/Q* 

which is still bounded in L°°(0, T; L P (U)) for all p < 00 . Therefore, in both cases 

He,5 


V$6 


is bounded in L°°(0, T; L p (fl)), p G [l,+oo). 


( 53 ) 


Control of the singular pressure. Since we have the relation 

KeA*s) = Hi, S (®S) 
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We can deduce directly controls on the singular pressure 


7 T £j , 5 (^* 5 ) is bounded in L°°(0, T; L P (Q )), p£[l,+oo). 


( 54 ) 


and 



bounded in L 2 (0, T; L q (Q)), q £ [1,2). 


(55) 


Remark: Compared to the work done with constant viscosities in [31], the singular viscosities 


via the K-entropy estimate (1341) provide directly an uniform control of the singular pressure 
without additional estimates using the Bogovskii operator. In addition, we have much more 
integrability in the present case thanks to the K-entropy which controls finally X7 p £: s(&) and 
consequently all the powers of fi e 5 ( 4 ?) and 7T e $($). In comparison, with constant viscosities we 
only get 7r £! <5(4>) bounded in L 1 ((0,T) x 12) which forces us to derive an additionnal estimate 
in order to prove the equi-integrability of the sequence. 

Convergences 

Following the classical steps of the stability of weak solutions of Navier-Stokes equations 
with degenerate viscosities (see for instance [28]). we prove first that the volume fraction 
4*5 converges to 4> £! $ a.e. and in C([0, T]; L P (Q)) for all p £ [l,+oo) thanks to the Aubin- 
Lions-Simon lemma (see [34]). This convergence leads then to the strong convergence of the 
pressures p( 4 >, 5 ) and ^^($, 5 ) in C([0, T], T 1 (fl)). 

Using again the Aubin-Lions-Simon lemma, we get the strong convergence in L 2 (0, T; L 9 (I2)) 
for all q £ [1,2) and the convergence a.e. of ms = Q 5 U 5 towards some m. Since ms/s/Qs is 
bounded in L°°(0, T; L 2 (U)), Fatou’s lemma provides 



which implies in particular that m = 0 a.e. on the vacuum set = 0}. Therefore we can 

define a limit velocity u e such that 


m(t, x ) 


if 

if 4^ = 0 


x) = < <S> Si &{t,x) 

0 


Finally &sus converges a.e. and strongly to in L 2 ((0, T) x 17). 

At the limit we recover the maximal volume fraction constraint 


Lemma 1 At the limit <5 —>■ 0, we have meas{(t, x) : $> £> tf(t,x) > 4>*} = 0. 
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Proof. This result is based on the control of the singular potentiel energy when 5 > 0 


T 5 e £j5 ($ 5 )l 

{<3><5—<5} 

' vr e , 6 (s) 

to 


> $«5 


ds 1 


(«1> 5 /<&* >1—<5} 


- (cj>*)2(l _ 5)2 \J Q 


$8 


r<5>*(l-5) 


7r e,5(s)ds l{$ 5 /$*>i-<5} 


> 


> 


> 


L {$W$>*>l-a} [ 
£<&*(1 - 5) 

£$*(1 — 5 ) y 


r$*(l-<5) 


10 


r>$*(l-<5) 


S \ 7 

<F 

S \7 

<F 


exp 


1 —< 3 . 


1 - s/<&* 


r lH ~CL 


1 - s/$* 


ds 


pi—5 

/ * 

Jo 


r 7+l 


dr 


1-5 J 0 1 - t 
> {Cl (— log 5) - C £ )l{$ s /$*>i-6} 

Integrating over 17 and letting <5 go to 0 we recover at the limit 

meas{(7, x) : <h £j ,?(t, x) > $*} = 0 . 



We can also prove the strong convergence of the singular terms p £ ^(d^), 7r e< j($,s) and 
^£, 5 (^ 5 ) in T p ((0,T) x 17) for all p £ [l,+oo) . Indeed we have previously seen that p £t s(^s), 
7r e:,<5(^< 5 ) and -^£, 5 (^ 5 ) are bounded in L p ((0, T ) x 17) for all p £ [1, + 00 ). Besides #5 converges 
a.e. and strongly to <h £i ,j and we ensure that meas{(f,x) : $> £j &(t,x) > <!>*} = 0. Therefore 
we guarantee that ^^($ 5 ), 7r £)5 (^ 5 ) and A e ^($ 5 ) converge a.e. towards p £ (® £ ^), 7r £ (<F £;1? ) 
and A £ (<f> £) 1 j) respectively. The Dominated Convergence Theorem finally provides the strong 
convergence of p £ ^, ir £t § and A £j g in L p (( 0, T) x 17) for all p £ [1, + 00 ). 


Convergence in the drag term. We want to show that 

r -T 


/ / |3\5KK “ -> 0 

Jo Jn 


For that purpose, we introduce R > 0 and split the previous integral into three parts 
f T 


/ / - $ £ , 0 |m £ , 0 |u £ , 0 | dxdt 

Jo Jn 

< / $ 5 |U( 5 |u 5 1 {| U5 |< R } 

Jo Jn 

+ / / |$5KKl{|u 4 |>fl}| dxdt 

Jo Jn 

+ n 


^ > £,^|' U '£,1?|^ t £,l?l{|u ej ^|<i?} 


dx dt 


/o Jn 

— rf» — 


£ ,l9 I U £, 1? I U £, 1? 1 {| U e: 0 I > R } 


dx dt 


First we have = <I > 5 which converges a.e. to ^ £ \^ 2 £ ^\u £ ^\u £ ^ on the set 

0} thanks to the convergence a.e. on (0,T) x 17 of &sus an d <3? < 5 . In addition 

|<I> ( 5 |ti 5 |u 5 l{| U5 |< R || < R 2 $s —> 0 on {<&£,,? = 0 } 
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Therefore we get the convergence a.e. of < F, 5 |« ( 5 |'U ( 5 to tf\u £ $\u e $ and the Dominated Con¬ 
vergence Theorem gives the convergence to 0 of the first integral, 


f f $5KKl{| Ui! |<.R} 

Jo Jn 


^£,1? | U E,fl \ U £,'& l{|« e ,,9 | <_R} 


dx dt —> 0. 


Concerning the two remaining integrals we use the control given by the K-entropy and we 
write 


^£,1? | ^£,1? | l{|u £ $ | >_R} 


dx d t 


/ / |$ 5 |« 5 |u 5 l{| U5 |> fl }| dxdt + 

Jo Jn Jo 

*1*51^.5| 3 1{|« 5 |>r} dx dt + / / ^ £) tf|u 6) tf| 3 l{| Uetf |>B}dxdi ) 

i Jo Jn J 


1 

< — 

“ R 

C 


f 

'0 Jn 


< 


R 


Letting R go to +oo, we obtain the strong convergence in L 1 ((0, T) x D) of the turbulent drag 
term towards $ £ ^\u ej #|u E y. 


Convergence in the convective term. For \fW 5 u 5 we develop the same idea as for the turbulent 
drag term and we decompose the integral between the small and the large velocities 


dx dt 


f [ \/$5US - V 7 $E,1?«£. 

Jo Jn 

< / / V / ^5 w <5 1 {|«d<- f ?} “ V®eJ> u e,V 1 {\ U£ #\<R} 

Jo Jn 

+ [ [ V^5US^-{\u s \>R} da 

Jo Jn 

+ / / y/^U £ ^l { \ u \> R} 

Jo Jn 


dx dt 


dx dt 
2 


dx dt 


As previously we can show the convergence a.e. of \fW$us to and the Dominated 

Convergence Theorem gives the convergence to 0 of the first integral. For the two last integrals 
we observe that we can split <f , |n| 2 as 

<f>|-u| 2 = <L 1/,3 $ 2 / 3 |n| 2 

and by the Holder inequality with 1 < q < 3/2 we get 

rT 


f [ $M 2 l{| Ui5 |xR}dxdf 
Jo Jn 


< 


< 


[ 

Jn 


/3 dx dt 


C 


R(3-2q)/q 


W / rT (■ \ 1 / c l 

j / <S> 2q/3 \u\ 2q dxdtj 
[ [ $|u| 3 dxdt) —► 0 

'0 Jn J R ^°° 


We conclude that 

V^sns converges strongly to y/^e ,$u e $ in L 2 ((0, T) x H). 
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Convergence in the diffusion terms. Since /u e 5 ($< 5 ) D (ug) and A e ^fH^div^ are bounded in 
L 2 (0, T; L 9 (12)) for all q € [1,2), we deduce that they converge weakly in L 2 (0, T; L 9 (12)) for 
all q € [1,2) towards /j, Ei g(^) D (u) and A £!< y($)divu respectively. 

Remember that the diffusion terms make sense in the weak formulation of the momentum 
equation if they are written under the form 


2 / [ he,s(® 8 ) D (ug) : VC 

Jo Jn 


n v** 

Jo Jn 

f [ y/** 

Jo Jn 


dihe,s(*&s') 

V®5 

djUeji®, 5 ) 

V&8 


diC j + 


di( j + 


he,s{^s) 

V$5 

he,5(^5) 


dk j 


dU j (56) 


For the first integral we have the strong convergence in L 2 ((0,T) X 11) of \fWfug towards 
and the weak convergence in L 2 ((0,T) x 12) of - towards ^ thanks 


V $8 




to the control (15311 . 

Recall that we define V £) g such that (1471) which allows us to write for the second integral 


= V$~8V<Pc,8(*8) = W e ,*($«) 

and we can prove that VI4 ^(S^) weakly converges in L 2 ((0, T ) x 12) towards W £ . By unique¬ 
ness of the limit in the sense of distribution we have y /< F e ,$V 14 = V// £ (3> £i tf). Thus we can 
pass to the limit in the weak formulation f[56[> to obtain (1221) . 

For the other diffusion term, 


[ f 7^== J M £ ,<5 (^ > <5)div n^div ( 
Jo Jn VheA^s) v 


we can prove the strong convergence in L 2 ((0,T) x 11) of the first part 


A £ ,<5( < 1 > <5) 

Ahe,s($s) 


Indeed 


we have the bound (1521) and the convergence a.e. since converges strongly to <h £) ^ and 
<h £;1? < 4>* a.e.. On the other hand fi e> s($>s)divug converges weakly in L 2 ((0,T) X 12) 
towards sj[i £ (<& ££ q)<\i\ u £ j) (this is the previous point). We deduce then the convergence of 
the integral towards 



Jhei^eA^ « £ ,^divC 


Note finally that, at the limit 5 = 0, we have the relation 


7Te ($£,#) 



7 

he (^e,i?)> 


(57) 


Remark on the limit passage 1 ? —> 0 : Since at this stage we ensure that $ £i ^ is bounded 
in L°°((0, T) x 11) we can deduce a control of V# £i ^ which does not depend on d. Indeed, 
thanks to the K-entropy inequality (1341) we have 

= V^V* e (<M bounded in L°°(0, T; L 2 (H)). 
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( 58 ) 


Then, since n' e (& e $) ^ 1 and € L°°((0,T) x P), we get that 

V$ e> # is bounded in L°°(0, T; L 2 (P)). 

We can then pass to the limit & —> 0 in the equations to obtain system dm 


4 Proof of Theorem [2], recovering the two-phase system as £ -A 0 


The aim of this section is to rigourously prove the limit passage from the suspension model 
m towards the two-phase system (12411 . What differs from the previous section is that at the 
limit volume fraction $ can reach the constraint <h* on a set of positive measure. We expect 
then that the bounds on the diffusion terms will be more subtle because the singular terms 
involve the quantity 

£ i+ a 



We see the competition between e 1+a which tends to 0 and 1 


(J>* 


which can tend to 0 possibly 


faster than e 1+a . In particular we do not have a uniform control on A e and to pass to the 
limit in the corresponding diffusion term we will need to consider the renormalized continuity 
equation. 


Estimates 


Control of the volume fraction. From the previous step, we know that is in L°°((0, T) x P) 
since it is bounded by <h*. As it has been explained for the limit passage i? —>• 0, we have in 
addition an uniform control of the gradient V<l? e in L°°(0, T; L 2 (P)) since 




< vT* 


\44M 

V $~6 




<E L°°(0,T;L 2 (P)). 


Integrahility given by the drag term. As previously we ensure extra-integrability of 
thanks to the the turbulent drag term present in the momentum equation 

$> x J^u e is bounded in L 3 ((0,T) x P). (59) 


Control of the singular coefficients 

vM$ e ) 

and 

J'Me(^e) 

since moreover <I> S < ( l>*. then 

Mc( < h £ ) ii 


Since the K-entropy m gives 
is bounded in L°°(0, T; L 2 (P)) 

is bounded in L oo (0,T;L 1 (P)), 

i bounded in L°°(0, T; W 1)2 (P)) 


28 






Thanks to the relationship between n e and /A we deduce that 


Vvr £ ($ e ) = 


7 ^r 1 

($*)7 


^(<l> e )V<& £ + ^ V^($ e , 




cj>=* 


is bounded in L°°(0,T; L q (Q)), V q € [1,2). 


( 60 ) 




As for the previous step <5 —>• 0, we need a control of the quantity f . Let <1> £ (0,1), 


on the set {<3? £ /<l>* > 1}, 


< Ve{$s) 


is bounded in L°°(0, T; L p (0)), p < oo. On the other set < <h}, <h £ is far from <L* and 

then uniformly in e 

-1-| -a 

exp 


/r £ ($ £ ) < vw; 


V®1 


1 - 


-1 )V®~e< Cy/¥ e 


which is still bounded in L°°(0, T; L p (fl)) for all p £ [l,oo). Therefore, in both cases 

/a($ £ 




is bounded in L°°(0, T; L p (fi)), p€ [l,+oo). 


(61) 


Passage to limit e —>■ 0 

With all these estimates we can now pass to the limit in the weak formulations of the mass 
and the momentum equations. The main convergence arguments remain the same as those 
presented in the previous section. We ensure with the Aubin-Lions-Simon lemma the strong 
convergence of the density <f> £ towards <f> in C([0, T], L p (£l)) for all p £ [1, +oo) and the limit 
gradient of the density Vd> is bounded in L°°(0, T; L 2 (fl)). We can also obtain the strong 
convergence of the momentum <L £ u £ towards m, define a limit velocity u equal to 0 on the set 
{$ = 0} and such that m = <!>«. 

The procedure to pass to the limit in the turbulent drag term and the convective term is 
the same as previously and consists essentially in splitting the integral (for instance of the 
convective term) 

I | y/W e u e ~ V$u\ 2 

Cl 

between small and large velocities. The Dominated Convergence Theorem gives the conver¬ 
gence to 0 for the small velocities whereas we need (15911 to control the integral for the large 
velocities. 

With the previous controls we get that 

' vr £ ($ £ )—>U weakly-* in L°°(0, T; W 1,I? (D)) 

< V 9 € [1, 2) 

k /A($ e ) = —> p + $ weakly-* in L°°(0, T; W 1 - 2 ^)) 
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Lemma 2 At the limit e = 0 we have 


and the equality 


In particular, 


($* - $)/Z = 0 
(<T - $)II = 0 



n e L°°(0, T ; VU 1,2 (fl)). 


(62) 

(63) 

(64) 


Proof. When e > 0 we have 


(l-$ £ /$>i($ £ ) = 


$e(l - 4> £ /$*) 


exp 


1 —|—CL 




‘ £ £*+“ 


1 - $ £ /$* 
* / / pl+a 

exp 


1 - 


- 1 

-11. 


We need to consider separetely three differents subsets 
S3 1 = {3b < a, s.t. V e, 1 — <3? e /4>* > s b }, fl 2 = {3c > 1 + a, s.t. V e, 1 — <l> £ /<f>* < e c }, 

n 3 = {V £, £ 1+a < 1 - $ £ /$* < £ a } 

• on fi 1 , this is the case where 


^exp 


T l-\-a 


-1 -> o 


i - $ £ /<i>* 

and for which we have directly the convergence of (1 — <3? £ /<f>*)/i{ to 0. 
on It 2 , the most singular case, we have 


1 —|—< 3 - 


X = 


+oo 


1 - <f> £ /<l>* 

and if p > 0 and X is large enough (or £ small enough), we can ensure that 

4 < (exp(X) - 1)». 

Using the fact that <1> £ is bounded away from 0, we get 

A+a \ \ 1+ P 


(1 - < £ a $s exp 


1 - 


- 1 


* 1+p 

< £ a + l +P($)-P £ 


£ £i+P 


exp 


1 —|—Q. 


1 - 3> £ /5>* 


- 1 


1+P 


= C£ a+1+p $- p (/r{(<h e )) 1+ P 

which tends to 0 since /r{(<3? £ ) is bounded in L°°(0, T; L p (Pl)) for all p £ [l,+oo). 
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• on II 3 , the intermediate case, we ensure that 


exp 


then, since we assumed that a > 1 




1 - < h £ /d> i| 


- 1 < C 


i 1 — <f> £ /<l>* 

(1 - * e / S >2 ^ C ^e Ja 


< Ce 


a—l 


o 


Thus, in every case, (1 — d> £ /<f>*)// £ converges strongly to 0 in L°°(0, T; L p (fl)) for all p E 
[1, Too). With the same arguments, since 7r e ($ £ ) = ^ /i £ (T £ ) has the same divergence 

close to $*, and since <h £ is in L°°(0, T; L p (fl)), for all p E [l,+oo), we ensure that 

(i - $/$*) n = 0 . 


By the strong convergence of <f> £ and the weak convergences of p £ (<& £ ) and tt £ ( < I ) £ ) we get in 
addition 



Combined with two previous constraints, II = p = 0 on {<h < < h*}, it gives hnally 


II = p, a.e. on (0, T) x Q 
and since Jl lies in L°°(0, T; VF 1,2 (fl)), 

n E L°°(0,T; W 1,2 (fi)). □ 


Finally 


tt £ ($ £ ) —> n weakly-* in L°°(0, T; W' 1 '«(fi)) 

V?G [1,2) 

. Pe($e) = nl($ e ) + $ £ ^ n + $ weakly-* in L°°(0, T; W 1 ’ 2 ^)) 


Concerning the diffusion term /x e ($ £ ) D (u £ ), the weak formulation writes as 

rT 


If 

rT 


JW £ u j (d 2 C j 

v e v 13 ^ ^ 


/ / V*e< 

Jo Jn 


i ( p £ {*& e) q ^-j _|_ djp £ (& £ ) g j 






The hrst integral converges to 











since 


M$e) 


converges weakly-* in L°°(0, T; L P (Q)) for all p G [1, +oo). We recall that 


djVs($ s ) 


djVe^e) 


converges weakly in L 2 ((0,T) x Q) and by uniqueness of the limit in the sense of distribution 

\/$W = v(n + $). 


Therefore the second integral converges to 

rT r ^(II + T) 


[ [ 

Jo Jn 


\/¥ 


di( j + 


fl(n + $) 

\/I> 



Transport equation relating II and A. 


Let us write the renormalized continuity equation on 


dtpli^e) + div (fil($ e )u E ) + 


M$ e ) 

2 


div u e = 0 


(65) 


or if we write the weak formulation 

rT 


f [ [ 

Jo Jn Jo 




\/W £ u £ ■ V£ 


In 

- ^(A £ (^ , e )div u £ , C) = J 0) 


for £ G D([0, T) X Q). Using the convergence already mentionned : /4(d> £ ) converges weakly-* 
in L°°(0, T; L P (Q)) for allp € [l,+oo) towards II, p\{<& £ )/\/WJ converges weakly in L 2 ((0,T)x 
U) towards II /\f$ and \/WJu £ converges strongly to \/l>u in L 2 ((0, T) x Q), we deduce that 
A £ ( ( L £ )divrt e converges in the sense of distributions towards a weak limit denoted by A and 
belonging to IU^ 1,oo (0, T; L p (I2))nL 2 (0, T; I / U _1,9 (H)) for all p G [1,+oo), <7 G [1,2). Moreover 
the equation (I24c() is satisfied in the sense of distributions 


d t U + div (Hu) 


A 

2 


This completes the proof of Theorem [2] 


5 Incompressible flows with pressure dependent viscosity 

This last section is devoted to the proof of Theorem [3] Our first goal is to show that the limit 
continuity equation (I24al) associated to the constraint 0 < < d>* is compatible with the 

incompressibility condition divu = 0 on the set {4? = 4>*}. Next we prove that the suspension 
model with initial density = $*(1 — £ a 4>*/II 0 ) approximates thanks to Theorem [2] the fully 
incompressible system with pressure dependent viscosity. 


We need to extend the compatibility lemma given by P.- L. LIONS and N. MASMOUDI in |24| 
to the case of degenerate viscosities. 
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Proposition 2 (Compatibility relation) Let (<h,u) be such that 

$eL p ((0,T)x!l) Mp € [l,+oo), V<I> € L°°(0, T ; L 2 (Ll)) 

\/$VueL 2 ((0,T)xn), $«£ L°°(0,T;L 9 (O)) with q > 1 

satisfying the continuity equation 

dt$ + div($u) = 0 m (0,T)xfl), $(0) = $°. 

TTien t/ie following assertions are equivalent 

1. divu = 0 a.e. on {4? > <!>*} and 0 < 4>° < 4>*. 

2 . 0 < 4> < <f>* 


Proof. 

• (1 =>- 2) As in [24], we set 


m 


0 if r < 0 

< r if 0 < r < <3?* 

1 if r > <f>* 


and Prj a regular approximation of f3 such that /3 v (r) = /3(r) on 

(—oo, —77) U (77, <E>* — 77) U ( 4 >* + 77, +00) 

and such that 

Since V<I> £ L°°(0, T; L 2 (fl)), we can multiply the continuity equation by /3'(<f>) and 
obtain 

dtPrii®) + div (/3 jj(3>u)) + ((£,,)'($)$ - /3 ?7 ($))divn = 0 ( 66 ) 

We have that /3 J? ( < f>) converges pointwise and in L 2 ((0,T) x fl) to /3(<f>). Moreover 

((/3 t? ) / (<P)<P — /^ (<!>)) div u = — \/^divn is bounded in L 2 ((0,T)xO) 

v ‘h 

and converges to l{<j>>$*}div 77 . Then passing to the limit in (16611 with respect to 77 and 
using the assumption divu = 0 on {<h > <3?*}, we get 

dt/3(&) + div (/3(4>)u) = 0. (67) 

To conclude we set d = /3(4>) — <!>, regularizing the function |d|, we show as previously 
that |d| satisfies 

f d t \d\ + div (\d\u) = 0 

l M(o) = 0 

Then integrating in space we get 

/ |d|(t)= [ |d|(0)=0 
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and therefore 

which means that /3(<1) = <f> or 


d(t) = 0 for all t 
0 < <f> < $*. 


(2 =>- 1) Assuming that 0 < < <3?*, equation (1671) holds for /3(4>) = 

integer k, since <h* is a constant 


<f> 

cj>* 


$ 


On the left-hand side we have 
dt 


$ 




d t( ^7 + div h^7 u = (! - k ) diYU 


$ 




, for any 


( 68 ) 


(|r) € , div ^ 

for all p 6 [1, oo) and for a q € (1, +oo), which shows that the right-hand side of (1681) is 
a bounded distribution. Then, if we let k go to +oo we obtain 

divu—tO in T>'({0 ,T) x fi). 

/ $ \ k 

On the other hand, I — j divrt converges pointwise to !{$=$*} diva and si 


since 


<l> 

—- ) div u 
$ 


fyk-l/2 


\/¥ |div it| is bounded in L 2 ((0,T) x Q) 


we conclude by uniqueness of the limit in the sense of distribution that 

l{ $=$ *}divrt = 0. □ 


5.1 Existence for the incompressible system with additional drag 

Let us prove now the Theorem [3] we consider for that the approximate initial data 


u° e eL 2 (n), <L° = <r 


(69) 


such that 


—> V&*u° in L 2 (fl) 


with e small enough to ensure 


1 -£ a - 


<£>* 


> 0 . 


ess inf II 0 

The approximate initial volume fraction obviously satisfies hypothesis (1151) . is positive thanks 
to the previous assumption and bounded uniformly with respect to e in IT 1,2 (11). One can 
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also check condition © , i.e. \m° £ \ since <f>0 is far from 0 and m £ = < 3?°k° is 

L 2 (Q). 

The approximate pressure 7r £ (<!>0) converges a.e. to II 0 > 0 since a > 1 and 


m 




< f >°\ 7 d >° 
d> 


1 —< 3 . 


= 


1 - <f>o/d* 

/ a <f*\ 7+1 exp( e n°/<h*) - 1 

V e W 


<f> 


- ( 1 - £“(7 + 1)|^ + o(e“) ) ( -j- + £7^2 + °( £ ) 


n° (n 1 


0\2 


$* (<J> 


,(n°) 


0\2 


= n °+ £A ^+ 0 ( £ ) 


Let us check now that m and f!18[) are also satished. We have 


VM<h°) = 


-I exp 


~l+a 


1 - 


- 1 + e a 


dL 






-I exp 


d>* (1 - d>°/d>*) . 

vn° 


£ n°\ \ $ e ( n °) 2 

<!>* J J <&* 


then we can bound 


|V/i e ($j?)l < 


$* 

rp 


(exp (n°/d>*) -1) +1 |vn°| 


Therefore, thanks to the control of n° in L°°(fi) HIT 1 ’ 2 (id) we deduce that // £ (d>0) is controlled 
in LF 1,2 (fd). Finally, since d> £ i s bounded by below, we check the condition (1171) 




\L 2 8 


< c 


Concerning condition (jl8|) . we establish the following control 

d>o e£ (d> £ °)=d>° r e ^-ds 

Jo 


d> 


0 /■$? 


*17 


e(d» 


$0 

< e 


*17 


e(d» 


$0 

< e 


e ($*)7 

'eU° 

exp 


L ' s7_1 ( exp (r 

r' s1_I ( exp (r 


^ 1 —)— 


- s/<&* 

crl+tt 


- *° e /$ 
eU° 


— 1 J ds 
— 1 | ds 


<t>* 


— 1 ds 


< 


d>* / 1 (d> 0) 7+1 


£7 


($*)7 


exp(n°/^)-l ^ 0 

7 £ 
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The last quantity is then bounded in L 1 ((0,T) x fi) as desired since II 0 is in L°°(Q) and 
is in L 1 (17). This proves the condition (1181) . 

If we consider now the solution (<3? e , u e , 7T e (3> £ )) of fllOal) f llObl) with initial data 
by the conservation of the mass we have 


1 - e a 




ess 


inf II 0 


l«l< 


$° = 




(70) 


and therefore 






(71) 


Jn 

Theorem |2] ensures the a.e. convergence of <I> £ towards a limit <h. Moreover, satisfies the 
constraint 


o < < <r 


(72) 


Then, necessarily by conditions (1711) ( 1721) we have $ = <h* a.e. on (0, T) X Q which proves 
with Proposition [2] the existence of a global weak solution (u, II, A) to 


div u = 0 


7 

Jo Jn 

r T 


(73a) 

(73b) 


mu- f f n u-vz+£,o= [ u°m 

[ [ u-dtC- [ [ (u®u) : VC~ 

J o Jn Jo Jn 

+(A idivC) + 2((n/$* + l)D(u),VC)= / «°-C(0), (73c) 

J 


I J i7 div ^ + r / / 
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